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We have studied the dynamical thermal stability of a one-dimensional topological insulator (TI) 
by deriving a master equation describing the coupling of the Tl-fermionic degrees of freedom to 
bosons in thermal baths at a given temperature. Unlike Hamiltonian perturbations or disorder 
effects, we find that the topological edge states are destroyed by thermal effects regardless the gap 
size of the Tl and for sufficient large times. Lifetimes for topological edge states are computed for 
physical implementations. Moreover, the non-equilibrium thermal evolution is able to distinguish 
topological and non-topological insulators. 
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1. Introduction. — The stability of topological phases of 
matter, also known as topological orders [I], against ther- 
mal noise has provided several surprising results in the 
context of topological codes used in topological quantum 
information [H[3]. However, very little is known about 
the behavior of a topological insulator (TI) subject to 
the disturbing thermal effect of its surrounding environ- 
ment. This is of great relevance if we want to address 
key questions such as the robustness of TIs to thermal 
noise, existence of thermalization processes, use of TIs as 
platforms for quantum computation, etc. Topological in- 
sulators have emerged as a new type of quantum phase of 
matter [H [5] that was predicted theoretically to exist [EJ- 
US] and has been discovered experimentally [T4HT6] . Ex- 
ploring the possible features and uses of TIs has become 
a very active interdisciplinary field. For this, knowledge 
about their stability under non-equilibrium thermal dy- 
namics is crucial in assessing the feasibility of proposals 
in quantum computation, spintronics, etc. 

In this work we present a first-principle calculation to 
test several thermal effects on a ID TI out of equilibrium. 
In order to achieve these goals, we need first to specify 
two choices: the type of TI and the type of thermal baths. 
As for TI, we work with the Creutz Ladder (CL) which is 
a paradigmatic example of quasi-one dimensional fermion 
system that exhibits the fundamental properties of TIs. 
Namely, localized states in the bulk of system and edge 
states in the form of zero-energy modes at the boundary 
[TT] ITS] . It is a crucial remark that the presence of these 
edge states is independent of whether the system size is 
finite or infinite. They constitute a clear signature of a 
TI in the case of the CL. Although the first experimental 
realizations of TIs are in 2D and 3D, the case of ID TIs 
also appears in the so-called 'Periodic Table' of TIs [T9l - 
121"] . Moreover, there are recent proposals to realize TIs 
in ID optical lattices [55], and in particular the CL [2"5] . 

As for the environmental quantum noise, we model it 
in the form of local bosonic thermal baths (see Fig. [IJ. 
This is rather natural, but novel since usually the bath 
and system degrees of freedom are taken to be of the same 
type. Here, we deal with a fermionic system but a bath 
of bosons. The reason for this choice is inspired by the 




FIG. 1: The TI system is pictured in (b) as a ladder of 
hopping spinless fermions with amplitudes K (horizontal and 
diagonal) and M (vertical). The fermions are coupled to a 
magnetic gauge field on the lattice that is perpendicular to 
the plaquettes. This is a lattice gauge theory known as the 
Creutz Ladder (CL) [171 I18j . The wavy lines at each sites 
indicate interaction with thermal bosonic baths. In (a), the 
interaction vertex of fermions with the bosonic bath is shown. 



traditional electron-phonon interaction in crystal solids 
|24j . The main difference is that our thermal baths are 
local in order to simplify its study. Moreover, this local- 
ity also fits into the traditional scheme of perturbing the 
global properties of a topological order by means of local 
external noise, as it is natural in topological quantum in- 
formation. If the CL is realized with optical lattices, then 
these local baths can be though of as external photons. 
Therefore, the meaning of the bosonic baths will depend 
on the specific realization we choose. 

A common belief in TI theory is that the gap defining 
the topological phase is enough to protect the system 
against interactions, disorder and even dynamical effects 
|25j . but the effects of dynamical thermal noise has not 
been addressed thus far: 

i/ We have shown that it does not hold for finite temper- 
ature effects: the TI order gets lost regardless of the gap 
size (see Fig. [2] and Fig. [I]). 

ii/ We have observed that the existence of topological or- 
der, strongly influences the system-bath interaction (Fig. 
[2]). In particular, for our ID model, the decoherence 
process remarkably depends whether the system is in a 
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topological phase or not. 

iii/ Notably, the interaction with bosonic thermal baths 
does not lead this system to the thermal state. How- 
ever the asymptotic state reached in the low temperature 
regime is close to it (see Fig. [3|. 

Our fundamental result is the derivation of the mas- 
ter equation (10), ( |TT| ) for a TI under a bosonic thermal 
bath from which all our results are derived. The total 
Hamiltonian of the problem considered reads as follows: 
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FIG. 2: Fidelities for the state of the CL with the initial 
Fermi Sea. Here we have taken 9 = n/2, T = IK and oj c = 
3/(2K) (of the order of the mean distance between the two 
bands). We see the fragility of the topological phase (m < 1) 
in comparison with the rest of the cases (m > 1). The main 
plot corresponds to the 8 sites CL, and the inset depicts the 
initial decay rate in the thermodynamic limit (see main text). 



H — H s + -ffbath + -Hint- 

The first term, H s , is the Hamiltonian of the CL, 
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here a n and b n are fermionic operators satisfying the anti- 
commutation relations: {a n ,a^/} = 5 n ,n', {&nj&„/} = 
5 n>n i. K and M are hopping amplitudes, while 9 is the 
magnetic flux per plaquette in natural units. It is known 
that for M < 2K and open boundary conditions, the sys- 
tem exhibits protected edge states [T7] that corresponds 
to a TI in ID fig]. 

The second term, -ffbath, is the free Hamiltonian of the 
local baths, 



(3) 



where A and stand for the bath bosonic operators that 
satisfy the canonical commutation relations [^,^t] = 



5 n ,n'fii,j , [A^jA^,] = 0. Moreover the index n denotes 
the position of the local bath on the CL, and i runs over 
the bath degrees of freedom. 

Finally, the third term in (JTJ), H lntl describes the in- 
teraction between the CL and the baths. In momentum 
space it reads (see Fig. [T]) 
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The quantity g l k k , regulates the boson-fcrmion coupling. 
Note that although the system contains free fermions in a 
gauge background field, its dynamics is highly non-trivial 
since the coupling with the bosonic bath involves three- 
body interactions (see Fig. [I]). 

2. Master Equation for a ID Topological Insulator. — The 
evolution of system and bath is given by the Liouville- 
von Neumann equation, which at the interaction picture 
reads (unless otherwise stated, natural units h = fee = 1 
are taken throughout the paper) 



where 
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Here f 1 := c and f 2 :— d are the operators which diag- 
onalize the CL Hamiltonian, i.e. H s = X^fceBZ ^i c \ c k + 
\2d\dk with 
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and m := M/2K. Moreover h% P k , := FgF%„ where 
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x k -(-l) a 
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Finally, := Ajj — are Bohr frequencies associated 
to the eigenvalues of the system Hamiltonian ^ which 
represent the two energy bands. 

By tracing out the bath's degrees of freedom from ([5| , 
we aim at writing a dynamical equation for the CL den- 
sity matrix, p s — Trb a th(p)- Under the natural assump- 
tions of the Born-Markov coupling to the thermal bath 
([53H25] and references therein), we arrive at the follow- 
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ing master equation for the TI: 
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These T's are the decay rates induced by the dissipa- 
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FIG. 3: Asymptotic occupation of the fermions for different 
values of m. The temperature is set to T = IK. On the 
left side we have represented some examples which present 
topological order, on the right side the system is out of the 
topological phase. Note that on the left side the point 4 is 
fixed whereas the point 5 is fixed on the right side. For the 
sake of comparison we have depicted also the Fermi-Dirac 
distribution which corresponds to the thermal state. 

tive dynamics in our system. Here, 0(w) denotes the 
Heaviside step function and n{ui) — [e w / T — l] is the 
number of bosons with frequency w in each local bath; 
T stands for the bath temperature. For the sake of sim- 
plicity we have assumed that g l k k , only depends on the 
difference between k and k' , g k k , = g\, where the energy 
e is related to k — k' through the dispersion relation of the 
baths, fn such a case, the so-called spectral density of the 
bath is formally written as J(ui) := J2i(9l) 2 ^( UJ ~ e% )- For 
dcfmiteness and as the CL could be realized in an optical 
lattice set up, we will consider a typical spectral density 
for a quantum optical ID system, J{uj) = acue^^ where 
a is a parameter that regulates the interaction strength 
(typically a will be the fine-structure constant) and uj c 
is a cutoff frequency. Furthermore, this "Ohmic" spec- 
tral density is widely used in the modeling of condensed 
matter systems as well [29] . 

Despite the apparent complicated structure of Y k ^ *| 
decay rates, the non-vanishing contributions are well un- 
derstood and both numerical and analytical calculations 
can be carried out with precision. 
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FIG. 4: Instability of the topological edge states against ther- 
mal noise in a 8 sites CL (the 6 sites CL provide same results) . 
Here 9 = n/2 and bath temperature is T = 0.04 K (i.e. 1% of 
the CL gap). The x— axis represent the longitudinal position 
on the CL. Thus, at t = all fermions are localized in left 
edge, and as time increases they tend to delocalize along the 
whole chain. The inset plot depicts the probability that the 
fermions are in the edge for several values of the magnetic 
flux 8. Thus it is clear that the topological edge states are 
unstable for any value of the magnetic flux. 



3.Non-Perturbative Thermal Dynamics. — In this section 
we analyze the physics out-of-equilibrium described by 



the master equation ( 10 ) in the case of a finite size CL. 
This allows us to retrieve non-trivial results about the 
stability of the topological order and whether it ther- 
malizes, among other properties. In order to study the 
stability of the system, we may chose different figures of 
merit. For instance, the fidelity T of the evolved mixed 
state of the system p s (t) and the initial Fermi Sea (FS) 
for the lower band of the TI, represents a measure of how 
the system remains correlated to its initial state which 
exhibits a topological order. 



7-[|FS),p 8 (t)] := (FS|e it£ (|FS)(FS|)|FS), 



(12) 



In fact, if this fidelity remains close to one then it is a 
strong indication that the topological order is preserved. 
Figure [2] shows the behavior of the fidelity in a CL of 
size N = 8. For some cases the fidelity may remain high, 
particularly for m > I, however in this case the CL is 
out of the topological phase [TT] [TS] . On the contrary, if 
m < 1 the fidelity may be reduced up to 10% of its initial 
value. 

A complementary criterium for the topological order to 
persist is given by the evolution of the fermion occupa- 
tion numbers. Most of fermions escaping from the lower 
band is a signature that the system no longer keeps the 
topological order. In figure[3]the occupation numbers are 
plotted for different values of m. The occupation is close 
to the Fermi-Dirac statistics, but they do not exactly fit 
each other. 

Finally, the existence of edge states is a well-defined 
property that characterizes a TI. If they disappear after 
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the TI is in contact with a thermal bath, we may unequiv- 
ocally conclude that the this type of topological order is 
lost. As shown in Fig. |4j they are unstable and tend to 
delocalize along the chain in time. Note that if the bath 
temperature is small in comparison with the gap of the 
CL, the system takes a lot of time to delocalize. This fits 
with the very well-known argument that TI insulators 
are stable to perturbations if they present a large gap. 
However, they always delocalize under thermal noise af- 
ter some sufficiently large period of time. 

Let us see the implications of our thermal evolution 
analysis on the physical implementations of the Creutz 
ladder with an optical lattice [23]. The thermal noise 
can be a model for heating induced by lasers that create 
the optical trap (due to fluctuating intensity profiles), 
or any other type of bosonic thermal noise. We take 
experimental values for to = 0, 9 = ir/2 as in [22/. K ~ 
3 h kHz, gap A = 12h kHz, and bath temperature T ~ 56 
nK which are currently reachable. We obtain a lifetime 
for edge states r ~ 67 ms . Hence, considering this type 
of thermal noise, measuring topologically ordered states 
could be possible within an optical lattice set up. 
^. Fidelity in the Thermodynamic Limit. — The evolution 
of the state can be written up to second order in time as: 
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Using ([12]), and after several calculations with the master 
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equation (10), we obtain the following result 
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At short times, two rates and T( 2 ) will determine 
how fast the fidelity of the Fermi Sea is lost during its 
evolution. The initial linear behavior for the lost of fi- 
delity given by is patent in Fig. [2] for N = 8 as 
well. Direct processes exciting electrons from one band 
to the other are dominant in the dissipative evolution 
as we might expect. Furthermore, in Fig. [2] inset, we 
can see that the initial decay of the Fermi Sea's fidelity 
strongly depends whether the system is in a topological 
phase m < 1 or not. More explicitly, the decay of fi- 
delity increases as we approach the topological crossover 
point m = 1, and then the decay decreases significantly 
for to > 1 - out of the topologically ordered regime -. 
This perturbative analysis for the thermodynamic limit 
is in total agreement with the exact results for TV = 8 as 
shown in Fig. [2] and with size N = 6 (not shown). 
5. Conclusions. We have derived a master equation 
describing the dynamical thermal effects of bosonic 
baths coupled to a one-dimensional TI. As this coupled 



fermionic-bosonic system is not exactly solvable, our for- 
malism is useful to address relevant thermal effects of 
TIs in ID. We observed that the dissipative dynamics 
distinguish whether the system is in a topological phase 
or not. We have shown that thermal noise delocalizes 
the topological edge states into the bulk bands of the 
TI for sufficiently large times and regardless of the gap 
size. While this is compatible with the existence of TIs 
in experiments, this thermal instability will play an im- 
portant role in detailed control manipulations needed for 
quantum computation. 
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